Black hole evaporation in a heat bath as a nonequilibrium process and its final fate 
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We consider a black hole in a heat bath, and the whole system which consists of the black hole 
and the heat bath is isolated from outside environments. When the black hole evaporates, the 
Hawking radiation causes an energy flow from the black hole to the heat bath. Therefore, since 
no energy flow arises in an equilibrium state, thermodynamic state of the whole system is not in 
an equilibrium. That is, in a region around the black hole, the matter field of Hawking radiation 
and that of heat bath should be in a nonequilibrium state due to the energy flow. Using a simple 
model which reflects the nonequilibrium nature of energy flow, we find the nonequilibrium effect on 
a black hole evaporation as follows: if the nonequilibrium region around black hole is not so large, 
the evaporation time scale of black hole in a heat bath becomes longer than that in an empty space 
(a situation without heat bath), because of the incoming energy flow from the heat bath to the black 
hole. However, if the nonequilibrium region around black hole is sufficiently large, the evaporation 
time scale in a heat bath becomes shorter than that in an empty space, because a nonequilibrium 
effect of temperature difference between the black hole and the heat bath appears as a strong energy 
extraction from the black hole by the heat bath. Further a specific nonequilibrium phenomenon is 
found: a quasi-equilibrium evaporation stage under the nonequilibrium effect proceeds abruptly to a 
quantum evaporation stage at semi-classical level (at black hole radius R g > Planck length) within a 
very short time scale with a strong burst of energy. (Contrary, when the nonequilibrium effect is not 
taken into account, a quasi-equilibrium stage proceeds smoothly to a quantum stage at R g < Planck 
length without so strong energy burst.) That is, the nonequilibrium effect of energy flow tends to 
make a black hole evaporation process more dynamical and to accelerate that process. Finally on 
the final fate of black hole evaporation, we find that, in order to make total entropy of the whole 
system increase along an evaporation process, a remnant should remain after the evaporation of 
black hole without respect to the size of nonequilibrium region around black hole. This implies that 
the information loss problem may disappear due to the nonequilibrium effect of energy flow. 
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I. INTRODUCTION 



Consider that a black hole is put in a heat bath, and that the whole system, which consists of the black hole 
and the heat bath, is isolated from outside environments. It has already been known that, by considering only the 
equilibrium states of the whole system,the equilibrium state is unstable for a sufficiently small black hole, and stable 
for a sufficiently large black hole [HQ. If the instability occurs for a small black hole and the system starts to 
evolve toward the other stable state, there are two possibilities of its evolution. The first is that, due to the statistical 
(and/or quantum) fluctuation, the temperature of heat bath exceeds that of black hole. Then the black hole swallows 
a part of heat bath and settles down to a stable equilibrium state of a larger black hole in a heat bath. The second 
possible evolution is that, due to the statistical (and/or quantum) fluctuation, the temperature of heat bath becomes 
lower than that of black hole. Then the black hole emits its mass energy by the Hawking radiation and settles down 
to some other stable state. However we do not know the details of end state of the second possibility, since the final 
fate of black hole evaporation is an unresolved issue at present. A more detailed explanation is given in section ITT1 

When one distinguishes the phase of the system by a criterion whether a black hole can exist stably in an equilibrium 
with a heat bath or not, the phase transition of the system occurs in varying the black hole radius. This phenomenon 
is known as the black hole phase transition [l[ [1| [I?} . This will be important in the study of primordial black holes 
Q, since the black holes would be formed in a radiation field in the early universe. And, in relation with the possibility 
of black hole creation in particle accelerators in the context of TeV gravity [ij , the black hole phase transition will be 
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important as well, since a black hole should be created in a quark-gluon plasma in our brane and in a graviton gas 
in the higher dimensional spacetimc. Therefore, although the black hole phase transition is an old and rather well 
studied topic, a more study on the black hole phase transition may give some contribution to the present trend in the 
study of black hole physics. 

As mentioned in the first paragraph, if an instability occurs for an equilibrium of a black hole with a heat bath, 
the black hole may evaporate. If the mass energy of black hole is radiated out completely by the Hawking radiation, 
then the thermal spectrum of Hawking radiation implies that only a matter field which is in a thermal equilibrium 
state may be left after the evaporation. This implies that the initial condition of black hole formation (gravitational 
collapse) is completely smeared out. For example, even if the initial state of collapsing matter is a pure quantum 
state, the final state after a black hole evaporation must be transformed to a thermal state. This contradicts the 
unitary invariance of quantum theory. This issue is known as the information loss problem || . The study on black 
hole phase transition, especially on a black hole evaporation after an instability of equilibrium occurs, may give some 
insight into the information loss problem. 

In this paper, we concentrate on a black hole evaporation in a heat bath after an instability of equilibrium occurs. 
During the evaporation process, there should exist an energy flow from the black hole to the heat bath due to the 
Hawking radiation. Since no energy flow arises in an equilibrium state, thermodynamic state of the whole system is 
not in an equilibrium. That is, in a region around a black hole, the matter field of Hawking radiation and that of 
heat bath should be in a nonequilibrium state due to the energy flow. However it should be noted here that the study 
on nonequilibrium phenomena is one of the most difficult subjects in physics, and it is impossible at present to treat 
the nonequilibrium nature of black hole evaporation in a full general relativistic framework. Hence we resort to a 
simplified model which reflects the nonequilibrium nature of black hole evaporation. Then we aim to find some insight 
into nonequilibrium effects on a black hole evaporation in a heat bath, and on the final fate of black hole evaporation. 

One may think that the generalised second law should be one of interesting issues concerning a black hole evapo- 
ration. For the case of a black hole evaporation in an empty space (a situation without heat bath), the generalised 
second law has already been discussed in reference @ with taking the nonequilibrium nature of black hole evaporation 
into account. The essence of the generalised second law in nonequilibrium situations has already been revealed in 
reference [6|. Hence the generalised second law is not in the scope of this paper. This paper may be considered as an 
extension of reference i 6] , although this paper is written independently of reference Q . 

The plan of this paper is as follows. Section UIl reviews the original discussion of black hole phase transition given 
in reference [l[ , which is based only on equilibrium states of the system. Then section IIIII is devoted to discussions 
about nonequilibrium effects of energy flow between a black hole and a heat bath, and insights into the nonequilibrium 
effects are presented. Summary and discussions are in section [TV] 

Throughout this paper, we use the Planck unit, h = c= G = ks = l- Then the Stefan-Boltzmann constant becomes 
a = 7r 2 /60, which is appropriate for a photon gas. When one consider a general non-self-interacting massless matter 
fields, as indicated in appendix IB! it is necessary to replace the Stefan-Boltzmann constant as a — ► a' = Ntt 2 /120 — 
N a/2, where N = + (7/8) n/, and rib is the number of helicities of massless bosonic fields and n/ is that of massless 
fermionic fields, (n& = 2 for photons.) 



II. EQUILIBRIUM MODEL 

Before proceeding to the study of nonequilibrium effects of energy flow between a black hole and a heat bath, 
we briefly review the original discussion of black hole phase transition given in reference This is based only on 
equilibrium states of the whole system which consists of a black hole and a heat bath. For simplicity we consider a 
Schwarzschild black hole in a heat bath. 

According to the black hole thermodynamics [3] [H , a Schwarzschild black hole can be modelled as a spherical black 
body whose equations of states are 



E Q = = — , S a = — L — = 7tR 2 a , (1) 

9 8irT g 2 ' 9 16ttT2 9 ' 



where R g , T gi S g are respectively the radius, temperature and entropy of the body, and E g is the internal energy 
of the body which corresponds to the mass of a black hole. It is obvious from these equations that the radius R g 
decreases when the body loses its energy E g . The black hole evaporation can be represented by the energy loss of 
this body. Further note that the heat capacity C g of this body is negative, 
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The negative heat capacity is a peculiar property of self-gravitating systems [9j]. That is, the energy E g includes 
self-gravitational effects of a black hole on its own thermodynamic state. Further it has already been revealed that 
fl(ij |. using the Euclidean path-integral method for a black hole spacetime and matter fields on it, the entropy of the 
whole gravitational field of a black hole spacetime is given by the quantity S g in equations of states ([T]). That is, the 
entropy S g includes not only the entropy of black hole but also the entropy of gravitational field outside the black 
hole. Hence it should be emphasised here that energetic and entropic properties of a black hole are encoded in the 
equations of states ([T]) , and we call the body the black hole. 

Put a black hole into a heat bath of heat capacity > 0, and set that the whole system which consists of the 

black hole and the heat bath is isolated from outside environments and in an equilibrium state ("micro-canonical 

ensemble"). The temperature of heat bath Th equals T g . We set C h = constant (> 0) for simplicity. The statistical 
(and/or quantum) fluctuation raises a temperature difference ST = T g — =/= 0. Then, there are two possibilities for 
the evolution of the system after the fluctuation occurs. The evolution is determined by the size of black hole and the 
signature of ST. 

Consider the case |C fl | > cf q) . If the initial fluctuation is ST > 0, the energy flows from the black hole to the 
heat bath due to the second law of ordinary thermodynamics, and relations dE g < and dEh > hold, where Eh is 
the energy of heat bath. Then dT g > and dTh > hold due to C g < and Ch > 0. Further, because the whole 

system is isolated, dE g + dEh = holds, and \dE g \ = dEh follows. Therefore we find dT g = dE g /C g < dEh/C^ = 
dTh dT g < dTh- This means that ST — > after a sufficiently long time. The same discussion holds for the initial 
fluctuation ST < 0. That is, the equilibrium state of the whole system is stable for a sufficiently large black hole 
which satisfies \C g \ = 2nR 2 g > C [ h eq) . 

Consider the case \C g \ < C ( h eq) . If the initial fluctuation is ST > 0, the energy flows from the black hole to the 
heat bath, dE g < and dEh > 0. Then dT g > and dTh > hold due to C g < and Ch > 0. However, the same 
calculations as in the previous paragraph result in inequality, dT g > dTh- This means ST — » oo and R g — > after a 
sufficiently long time, which corresponds to a black hole evaporation. On the other hand, for the initial fluctuation 
ST < 0, we obtain dT g < 0, dTh < and \dT g \ > \dTh\- This means that the black hole grows dR g > 0, and finally 

the heat capacity becomes to satisfy the stable condition \C g \ > C^ q \ In summary for the case \C g \ < C^ q \ the 
whole system transfers to a stable equilibrium state with sufficiently large black hole for a fluctuation ST < 0, or to 
some other stable state without a black hole (end state of black hole evaporation) for a fluctuation ST > 0. That is, 
the equilibrium state of the whole system is unstable for a sufficiently small black hole which satisfies | C' g \ < C h eq ^ . 

Hence from the above, we find that the phase transition of the whole system occurs in varying the black hole radius, 
where the phase of the system is distinguished by the criterion whether a black hole can exist stably in an equilibrium 
with a heat bath or not. This phenomenon is known as the black hole phase transition. 

The above discussion is primitive. Although a self-gravitational effect of a black hole on its own thermodynamic 
state is encoded in equations of states (P) , the above discussion ignores gravitational interactions between a heat bath 
and a black hole. However, using the Euclidean path-integral method for a black hole spacetime and matter fields on 
it, reference Q has extended this primitive discussion to include the gravitational interactions, then the same result 
(stable equilibrium for a large black hole and unstable equilibrium for a small black hole) has been obtained, although 
reference [2| considers the case that the temperature of heat bath is kept constant ( "canonical ensemble" ) . Here note 
that thermodynamic properties of an equilibrium state can be derived by statistical mechanics with using either a 
micro-canonical ensemble or a canonical ensemble. Hence, the difference of ensembles does not become an essential 
difference between references Q and @. That is, the primitive discussion given in reference [I[ retains the essence 
of black hole phase transition. An intuitive understanding why reference [l[ retains the essence is given in appendix 
|A"1 We can find an understanding by references [l[ and 0] and by appendix it is not the gravitational interactions 
between black hole and its environments but the self-gravitational effects of black hole on its own thermodynamic 
state (e.g. the negative heat capacity) that causes the black hole phase transition. In next section, we put our basis 
on reference [l[ in taking the nonequilibrium nature of black hole evaporation into account, and consider an isolated 
case of the whole system. 

III. NONEQUILIBRIUM EFFECTS 

A. Simple model of a black hole in a heat bath 

The discussion so far has been based only on equilibrium states of the whole system which consists of a black hole 
and a heat bath. However as mentioned in section |TJ there should be a nonequilibrium region around a black hole 
due to the energy flow which is caused by the energy exchange between the black hole and the heat bath. The energy 
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FIG. 1: NPT model. Left panel shows the particle flux between a black hole and a heat bath. Right one shows some variables 
and particles of radiation fields emitted by black hole and heat bath. 



exchange is never carried by a heat conduction but carried by an exchange of particle fluxes, since a black hole is 
not composed of matters but is a certain spacetime region [111 ]. The particle flux from the black hole is the Hawking 
radiation, and that into the black hole is the accretion of matters emitted by the heat bath. Hereafter, for simplicity 
and as explained at the end of appendix \X\ we consider general non-self-interacting massless matter fields as the 
matters which carry the energy exchange between the black hole and the heat bath. For example, photon, graviton, 
neutrino (if it is massless) and free Klein-Gordon field (□$ = 0) are candidates of such matter fields. Note that these 
fields possess the generalised Stefan-Boltzmann constant a' as mentioned at the end of section [JJ 

It is very difficult at present to treat nonequilibrium phenomena around a black hole in a full general relativistic 
framework. Therefore, in order to pick up the nonequilibrium effect of energy exchange between a black hole and a 
heat bath, we use a simple model which is obtained by modifying the equilibrium model considered in section [Til We 
remove a shell region around the black hole from the heat bath. Then the space of the removed region between black 
hole and heat bath is filled up with matter fields emitted by the black hole and the heat bath. Here note that, as 
explained in appendix when we use thermodynamic quantities which are evaluated at asymptotically flat region 
on a black hole spacetime, the thermodynamic property of the whole system which consists of a black hole, a heat 
bath and matter fields between them may be described by the model which ignores gravitational interactions among 
a black hole, a heat bath and matter fields. Hence, we consider the model named NPT after the Nonequilibrium black 
hole Phase Transition. 

NPT model: Put a spherical black body of temperature T g in a heat bath of temperature Th{^ T g ), and set that the 
whole system is isolated from outside environments. Consider the case that the equations of states of the black 
body are given by equations |T]). and that the black body and the heat bath emit some non-self- interacting 
massless matter fields with Planckian distribution (we call these fields the radiation fields hereafter). Then 
hollow a spherical shell region out of the heat bath around the black body as shown in figure [JJ The outermost 
radius Rh of hollow region satisfies Rh > R g by definition and the hollow region is concentric with the black 
body. The hollow region is filled up with the radiation fields emitted by the black body and the heat bath. 
Further we put two assumptions as follows. 

Fast propagation assumption: The volume of hollow region is not so large that the speed of light is approximated 
as infinity. Then the retarded effect on radiation fields during propagating in the hollow region is ignored. 

Quasi-equilibrium assumption: The time evolution in the NPT model is not so fast that thermodynamic states 
of black body and heat bath at each moment of their evolution are approximated well by equilibrium states 
individually. (Recall a quasi-static process in ordinary thermodynamics.) This is consistent with equations of 
states (JTJl which is obtained from a static black hole solution of the Einstein equation (Schwarzschild black hole). 
Further, since a Schwarzschild black hole is not a quantum one, the order relation R g > 0(1) should be required 
as well. 

Here note that, as mentioned in equation @, the energetic and entropic properties of a Schwarzschild black hole are 
encoded in equations of states (JTJ). Henceforth, as done in reference we call the black body the black hole. The 
reason why we consider an isolated case of the whole system is that we put our basis on the model used in section UTI 
and reference as mentioned at the end of section UTI 

There are three points which have to be mentioned here. The first point is about the thermodynamic states of 
black hole and heat bath under the quasi-equilibrium assumption. This denotes that time evolution of black hole is 
expressed as that thermodynamic state of black hole changes along a sequence of equilibrium states in the state space 
during time evolution of the whole system. And the same is true of time evolution of thermodynamic state of heat 
bath. Then the temperatures T g and Th are given as equilibrium temperatures. Therefore it is reasonable to consider 
that, during the time evolution, the temperatures T g and Th are constant within a time scale that one particle of 
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radiation fields travels in the hollow region until absorbed by black hole or heat bath. This is consistent with the fast 
propagation assumption. 

The second point is about the thermodynamic state of radiation fields under the quasi-equilibrium assumption. In 
the hollow region, the radiation fields of different temperatures T g and Th are simply superposed, since the radiation 
fields are non-self-interacting (a gas of collisionlcss particles). That is, the radiation fields are in a two-temperature 
nonequilibrium state. Here note that, since the temperatures T g and Th are constant within the time scale that one 
particle of radiation fields travels in the hollow region, it is reasonable to consider that the radiation fields have a 
stationary energy flow from the black hole to the heat bath (or from the heat bath to the black hole) within that time 
scale. Hence, at each moment of time evolution of the whole system, the thermodynamic state of radiation fields is 
well approximated to be a macroscopically stationary nonequilibrium state, which we call the steady state hereafter. 
The time evolution of radiation fields is expressed as that the thermodynamic state of radiation fields changes along 
a sequence of steady states in the state space. Therefore we need a thermodynamic formalism of two-temperature 
steady states for radiation fields in order to analyse the black hole phase transition and evaporation in the framework 
of NPT model. The steady state thermodynamics for radiation fields has already been formulated in reference [13] , 
and its brief summary is in appendix |B] We will make use of the steady state thermodynamics in analysing the NPT 
model. 

The third point which should be mentioned here is about the value of N which appears in the generalised Stefan- 
Boltzmann constant a' = Nir 2 /120. As mentioned at the end of section |T] and appendix [Bl N denotes the number of 
independent states in radiation fields. Further, because of the quasi-equilibrium assumption, the NPT model should 
be valid at most until the semi-classical stage of black hole evaporation ends (T g < Planck temperature) and no exotic 
particle is created in the radiation fields by the Hawking effect. Therefore it is reasonable to estimate the value of N 
by the number of independent states of standard particles (three generations of quarks and leptons and their internal 
states, and gauge particles of four fundamental interactions and their helicities). That is, we require the following 
condition hereafter, 

N = 0(10) (=► a' =0(1)) . (3) 



B. Energy transport and the phase transition in NPT model 

We discuss energetics of NPT model. The total energy of the whole system is E to t — E g + Eh + E ra d, where 
Eg is the energy of black hole given in equations of states |T]), Eh is the energy of heat bath defined by ordinary 
thermodynamics, and E ra d is the steady state energy of radiation fields in the hollow region. The energy E ra d is given 
by reference [l2j or equations (|B5j) and (|B8|) in appendix [5] (the generalised Stefan-Boltzmann constant a' appears 
in the form of E ra d, since the radiation fields are non-self-interacting massless matter fields), 

E r ad — 4(t' (G g + Gh 7ft) , G g = [ dx 3 g g (x) , G h = f dx 3 g h {x) , (4) 

Jv rad Jv rad 

where g g {x) is the solid-angle (divided by An) covered by directions of particles which are emitted by the black hole 
and come to a point x (see figure [1] and figure [5] in appendix [B)l . and gh(x) is defined similarly by particles emitted 
by the heat bath. Note that g g (x) + gh{x) = 1 holds by definition, and consequently V ra d — G g + Gh is the volume 
of hollow region. Further, since the black hole is concentric with the hollow region, we obtain g g (x) = ( 1 — cos8 r ) /2 
and gh(x) = ( 1 + cos6> r ) /2, where 9 r is the zenith angle which covers the black hole at a point of radial distance r 
(see right panel in figure [I}. Then the quantities G g and Gh are calculated, 
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Rl-R 3 g + (Rl-Rlf 2 } . (5) 



In order to understand the energy flow in NPT model, we divide the whole system into two sub-systems X and 
Y as follows. The sub-system X is composed of a black hole and "out-going" radiation fields emitted by the black 
hole, and the sub-system Y is composed of a heat bath and "in-going" radiation fields emitted by the heat bath (see 
left panel in figure [T]). That is, the sub-system X is a combined system of components of NPT model which share 
the temperature T g , and Y is that which share the temperature Th- Then the total energy of the whole system is 
expressed as 

Etot = Ex + Ey , (6) 
where Ex and Ey are respectively the energies of sub-systems X and Y, 

[E x =E g +E^ d , E^ d = Aa'G g T± 

\E Y = E h +E^ d , E^ d = Aa'G h Tf: ' 1 > 
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where E^ d and E^"' d are respectively the energies of out-going and in-going radiation fields. The energy flow in 
NPT model can be understood as an energy transport between sub-systems X and Y. In order to express the energy 
transport explicitly, recall that the energy exchange between a back hole and a heat bath is carried by the exchange 
of radiation fields emitted by them, as mentioned at the beginning of section fill Al Hence the Stcfan-Boltzmann law 
works well as an explicit expression of energy transport (with the generalised Stefan-Boltzmann constant a'), 
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where A g — Ai:R 2 g is the surface area of black hole, and t is a time coordinate which corresponds to a proper time 
of a rest observer at asymptotically flat region if we can extend the NPT model to a full general relativistic model. 
Note that, some of particles of the radiation fields emitted by heat bath are not absorbed by the black hole but return 
to the heat bath (see right panel in figure [1]) . Therefore the effective surface area through which the sub-system Y 
exchanges energy with the sub-system X is equal to the surface area of black hole A g , and the area A g appears in 
both of equations (jSJ . Further we note that the energy transport (jHJ) is formulated to be consistent with the isolated 
setting of the NPT model (E to t = constant) . 

It is useful to rewrite the energy transport ([8]) to a more convenient form for later discussions. By equations [[2]). 
(0) and jj}, the energy transport ([8]) becomes 
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where C g = —2nRg is given by equation j2]) and it is assumed for simplicity that Rh = constant and that Eh depends 
on Th but not on T g . The quantity Ch is the heat capacity of heat bath, and we assume Ch = constant > for 

simplicity. The quantity Cy is the heat capacity of sub-system Y under the change of T g with fixing T h , and Cy ^ 

is that under the change of Th with fixing T g . The quantity C^ d is the heat capacity of out-going radiation fields, 
and Cx is the heat capacity of sub-system X. Here note that Ex has no T^-dependence. In analysing the nonlinear 
differential equations ©, behaviours of various heat capacities (fTTj) are used. Some useful properties of these heat 
capacities are explained in appendices O El El an d El 

Here it should be pointed out that an inequality Cx + Cy^ < has to hold in order to guarantee the validity of NPT 
model. To understand this requirement, consider the case T g > Th for the first. The energy flows from a black hole 
to a heat bath via radiation fields, dE g < and dEh > 0. Then dT g > and dTh > hold due to C a < and Ch > 0. 
Recall that we are considering an isolated case of the whole system, E to t = constant 
Therefore, because of Cy^ > by definition, it is concluded that the inequality Cx 

inequality Cx < follows immediately due to Cy ' > by definition. The similar discussion holds for the case 
T g < Th, and gives the same inequality. Therefore the following inequality must hold in the framework of NPT model, 



(Cx+C { y 9) )dTg+C y y> dT h = 0. 
■f C^ < must hold. And an 



Cx + C^KQ (^C x <0). 



(12) 



This inequality is the condition which guarantees the validity of NPT model. A more detailed property of the combined 
heat capacity Cx + Cy^ is explained in appendix El which shows that inequality (TT2|) can hold for a sufficiently small 
Th- Therefore we assume that Th is small enough so that this validity condition (fT2"|) holds. 

Concerning the validity of NPT model, it is also important to consider what the quasi-equilibrium assumption 
implies. This assumption requires that the time evolution is not so fast. Therefore it should be required that, cither 
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when a black hole evaporates (T g > Th) or when it grows (T g < Th), the shrinkage (or expansion) speed of the surface 
of black hole is less than unity, \dR g /dt\ < 1. Using equations this becomes 



dR g 



dt 



1 

< 1 T 



T*C X 



<1. (13) 



This inequality is also the condition which guarantees the validity of NPT model. Hence, in the framework of NPT 
model, the analysis should be restricted within the situations which satisfy the conditions (JTSJ) and Q13p . 

Here it is helpful for later discussions to consider what a violation of validity conditions (fl2|) and (fl~3|) denotes. 
Firstly consider if the validity condition (fT2|) is not satisfied. Then the system, especially the radiation fields, can 
never be described with steady state thermodynamics ("steady" means "stationary nonequilibrium" ) . That is, the 
radiation fields should be in a highly nonequilibrium dynamical state. Then the quasi-equilibrium assumption is 
violated, because it is this assumption that lead us to utilise the steady state thermodynamics as explained in section 
MI Al Therefore highly nonequilibrium radiation fields make a black hole dynamical, and the black hole can not be 
treated by a stationary solution of the Einstein equation. Next consider if the validity condition (|13|) is not satisfied. 
Then a black hole evolves so fast that the quasi-equilibrium assumption is violated and that the black hole can not be 
described by a stationary solution of the Einstein equation. That is, a black hole should be described as a dynamical 
one, and consequently radiation fields evolve into a highly nonequilibrium dynamical state. Hence it is expected that, 
when one of the validity conditions Q12[) or (|13[) is violated, the system evolves into a highly nonequilibrium dynamical 
state which can not be treated in the framework of NPT model. 

Finally in this section, we discuss about the black hole phase transition in the NPT model. If the whole system is 
in an exact equilibrium, then a black hole, a heat bath and radiation fields have the same equilibrium temperature 
T eq (— T g = Th), and, as shown in section |TT1 one may think that an instability of equilibrium occurs when an 

inequality \C g \ < C^ e<? ' holds, where C h eq ^ = Ch + l6a'Tg q V ra ,i and 16 a' Tjf 9 V ra d is the heat capacity of radiation 
fields in an equilibrium. However we should be more careful to consider the criterion of occurrence of instability of 
the whole equilibrium state in the framework of NPT model. At the occurrence of instability, the temperature of 
black hole T g and that of heat bath Th become different from each other due to the temperature fluctuation. Then, 

referring to the conservation of total energy, dE to t — {Cx + Cy ) dT g + Cy ^ dTh — 0, a relation dT g > dTh holds 
when the following inequality is satisfied with the fluctuated temperatures T g and Th (^ T g ), 



C x + 4 9) 



<4 ft) => \C g \<cl 9 J d + Cy+C Y h \ (14) 



where > is used which is shown in appendix [E] This inequality denotes a growth of temperature difference 
ST = T g — Th- That is, we can expect that the criterion of occurrence of instability is given by inequality (| 14|) . 
However it should be noted that the temperatures T g and Th evolve in time for the case T g ^ Th- Therefore, even if 
inequality (|14p is not satisfied with specific values of fluctuated temperatures T g and Th at an initial time, there is 
a possibility that this inequality comes to be satisfied and an instability occurs during time evolution of the system 
(an example is shown in next section). And a converse situation is also possible that this inequality comes not to be 
satisfied even if it is satisfied at an initial time. Hence, when only initial values of fluctuated temperatures T g and Th 
are used, we can not exactly judge by inequality (|14[) whether an instability of the whole equilibrium state occurs or 
not. That is, in a practical usage, inequality (fT4"|) is not useful as a criterion of occurrence of instability. Here it should 
be noted as well that, contrary to the NPT model, the criterion of occurrence of instability in the equilibrium model 
used in section[H]is given by the inequality \C g \ < C^ 9 ' evaluated at the initial time. Hence the nonequilibrium effect 
of energy exchange tends to make the situation more complicated. 

C. Black hole evaporation in a heat bath 

1. General aspect of the NPT model 

In order to discuss a black hole evaporation after the instability of an equilibrium of the whole system occurs, we 
assume T g > Th hereafter. In order to analyse energy transport equations ^ from energetic viewpoint, we consider 
the energy emission rate J np t by black hole in the framework of NPT model, 

J npt = -^=a>§L(T*- T £)A g , (15) 
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where equations ([2]) and ([9]) are used. The larger the value of J np t, the more rapidly the energy E g of black hole 
decreases along its evaporation process. In other words, the stronger emission rate J npt denotes the acceleration of 
black hole evaporation. 

Here recall that, as mentioned in equations (jlip . we assume Rh — constant for simplicity as the evaporation 
process proceeds. However, because Rh is the parameter which controls the size of nonequilibrium region around 
black hole, it is useful to compare two situations which differ only by the value of Rh with sharing the same values 
of the other parameters of NPT model, R g , Th, Ch and N in a' . In order to make this comparison, we note the 
following three points; firstly C g < by definition, secondly Cx < shown by inequality (fT!?]) , and finally that \Cx\ 
is monotone decreasing as a function of Rh for Rh > (3V2/4) R g ~ 1.06R g (see appendix [D]) . The first and second 
points denote C g /Cx > 0, then the third point concludes that C g /Cx is monotone increasing as a function of Rh for 
Rh > (3-\/2/4) R g . Hence it is recognised that, for the case of Rh > (3\/2/4) R g , the larger we set the nonequilibrium 
region, the faster the black hole evaporation process proceeds. Numerical examples are shown later in figure [H 

The above discussion is a comparison of NPT model of a certain value of Rh with that of a different value of 
Rh. In the following subsections, we compare the NPT model with the other models of black hole evaporation, the 
equilibrium model used in section [TT] and the black hole evaporation in an empty space (a situation without heat 
bath). 



2. Comparison with the equilibrium model used in sectionW\ 

Since the NPT model is constructed from the equilibrium model used in section [Til it is interesting to compare the 
NPT model with the equilibrium model. The energy emission rate J eq by black hole in the equilibrium model is given 
with setting Rh — R g in J np t for all time, 

Jeq — 0' (T* — Tfl) A g = J , (16) 

where J is given by equation (flO|) . Then we find J npt = (C g /Cx) Jeq- Note that appendix [El shows C g /Cx > 1. 
Therefore, when the values of R g , Th and N in a' are the same for the NPT and equilibrium models, then J np t > Jeq 
holds. This implies that the black hole evaporation in NPT model proceeds faster than that in the equilibrium model. 
That is, we can recognise that the nonequilibrium effect of energy exchange between a black hole and a heat bath 
accelerates the black hole evaporation. 



3. Comparison with the black hole evaporation in an empty space 



Usually in many papers on black hole physics, the evaporation time scale of black hole is estimated with assuming 
that a black hole is in an empty space (a situation without heat bath). It is useful to compare the NPT model with 
the black hole evaporation in an empty space. When gravitational interactions between a black hole and matter fields 
of Hawking radiation are ignored (see appendix [XJ, the energy emission rate J empty by black hole in an empty space 
is given by simple Stefan-Boltzmann law in an empty space (dE g /dt = —cr'TgA g ), 

Jempty — & Tg Ag , (17) 

where it is assumed that matter fields of Hawking radiation is non-self-interacting massless matter fields (this is the 
same "radiation fields" as those considered in the NPT model). Then we find 

Jnpt = ~ J^j '^ em V t V ■ (18) 

Recall that T g > Th 1 — T^/T g < 1) holds generally for a black hole evaporation, and C g /Cx > 1 holds in the 
framework of NPT model (see appendix[E|. Then the factor (C g /C x ) (l - Ttl T f) 

may be greater or less than unity. 
Therefore it is not definitely clear which of J np t and J em pty is larger than the other. 

One may naively expect that the incoming energy flow from heat bath to black hole in NPT model never enhance 
the energy emission rate by black hole, and that the relation J npt > Jempty is impossible but J np t < Jempty must 
hold always. It is always true if the out-going energy flow due to Hawking radiation is never absorbed by the heat 
bath. However in the NPT model, the energy emitted by Hawking radiation is absorbed by the heat bath and affects 
the incoming energy flow from heat bath to black hole. Then it is reasonable to expect that the energetic interaction 
(energy exchange) between black hole and heat bath determines the energy emission rate J np t- When we take the 
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energetic interaction into account, a relation J npt > J empty which is naively unexpected is also possible as discussed 
in the following paragraphs. 

In order to discuss the energy emission rate J np t under the effects of energy exchange between a black hole and a 
heat bath, it is useful to recall the decomposition of the whole system of NPT model into the sub-systems X and Y, as 
considered in section UlI Bl Further it should be noted here that, from energetic viewpoint, a black hole evaporation in 
an empty space can be thought of as a relaxation process of an isolated system in which the sub-system Y is removed 
from the NPT model and the sub-system X is isolated. Therefore, for the black hole evaporation in an empty space, 
the energy emission rate J em pty is the energy transport just inside the sub-system X (from black hole to out-going 
radiation fields), and no energy flows out of the sub-system X. However the energy transport ([8]) of NPT model is 
the energy exchange between sub-systems X and Y. That is, in the NPT model, the energy Ex of sub-system X 
is extracted by the sub-system Y due to the temperature difference T g > T h , and energy flows from X to Y. The 
black hole evaporation in NPT model is essentially different from the black hole evaporation in an empty space due 
to the energetic interaction (energy exchange) between sub-systems X and Y. This difference may be understood 
by considering a limit of the energy transport as follows: one may expect that the energy emission equation of 
black hole evaporation in an empty space, dE g /dt = —Jempty, should be obtained from equations ([8]) by the limit 
operations, Th — > 0, Eh — > (remove the sub-system Y) and Rh — * oo (infinitely large volume of out-going radiation 
fields). However these operations transform equations {SJ into the set of equations, dE g /dt — — Jempty and = Jempty 
This contradicts the "evaporation", since an unphysical result E g = constant (= oo) is obtained. This implies that 
the black hole evaporation in an empty space can not be described as some limit situation of the NPT model. 

Hence, in addition to the naive expectation J npt < Jempty , the opposite expectation J npt > Jempty may be expected 
due to the energetic interaction as follows: because the energy source of out-going radiation fields is the mass energy 

Eg of black hole, it seems that, the more amount of energy E^L of out-going radiation fields is extracted by the 
sub-system Y, the more amount of mass energy E g of black hole should be radiated by the Hawking radiation. That 
is, when the energy Ex of sub-system X is extracted by the sub-system Y during a black hole evaporation process 
in NPT model, then the energy emission process by black hole (the black hole evaporation) is accelerated due to the 
energy extraction by sub-system Y. This implies J np t > Jempty 

The above discussion can be supported by the following rough analysis. When a black hole evaporates in NPT 
model, the temperature difference ST = T g — Th should grow infinitely, ST — > oo (recall the discussion in inequality 

(TH)) ). Then, because of equation (fT5|) together with the facts 1 — {Th/T g ) A — > 1 (as ST — > oo) and C g /Cx > 1 (see 
appendix [Ej, the larger the temperature difference ST, the larger the ratio J n pt/ Jempty Hence for the case of black 
hole evaporation in the NPT model, it is expected that the relation J npt > Jempty comes to be satisfied during the 
evaporation process even if the relation J npt < Jempty holds at initial time. And, if this relation J npt > Jempty holds 
for a sufficiently long time during the evaporation process, the evaporation time scale in NPT model can be shorter 
than that in an empty space. That is, it is possible that the black hole evaporation in NPT model proceeds faster than 
that in an empty space, where black holes of the same initial radius are considered in both cases. In next subsection, 
numerical examples which support this discussion are shown. 

4- Numerical example 

Since an analysis of nonlinear differential equations is difficult, it is helpful to show numerical solutions T g (t) and 
Th{t) of energy transport equations ©. The initial conditions are 

R g (0) = 100 (=S> T g (0) ~ 0.00079) , T h (0) = 0.0001 , (19) 

and the other parameters are set as 

C h = 1000 , N = 10 , (20) 

where see equation ([3]) for the value of N. Further we have to specify the outermost radius Rh of hollow region. As 
mentioned in equation (|15p . by the comparison of a numerical solution of energy transport ([9]) of a certain value of 
Rh with that of a different value of Rh, we can observe the nonequilibrium effect of energy exchange between a black 
hole and a heat bath. The numerical results are shown in figure [H and the value of Rh is attached in each panel. 
Time coordinate r in this figure is a time normalised as 



(21) 
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FIG. 2: Numerical solutions of energy transport (|9j for Ch = 
Th(0) = 0.0001. Horizontal line denotes the normalised time r = 
the outermost radius of hollow region Rh . 



10 , N = 10 and the initial conditions, R g (0) = 100 and 
t/tempty The size of nonequilibrium region is controlled by 



where t empty is the evaporation time (life time) of a black hole in an empty space (a situation without heat bath), 



which is determined as 
dE„ 



dt 



-J, 



empty 



Rg(t)=R g (0) l 
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(22) 



where one of equations U) is used, and t emp t y — 4.02 x 10 s by conditions (fT9"j) and (j2"D|) . This t empt y is usually adopted 
as the time scale of black hole evaporation in many papers on black hole physics. 

Further hereafter, we introduce the other time scale t npt which denotes a time at which one of the validity conditions 
of NPT model (JUJ) or is broken, 



L-npt 



t^ , to 



C x {ti) + C^\t l ) = Q,v{t 2 ) = 1 



(23) 



where Cx + Cy^ is treated as a function of time t through R g (t), and v(t) = \dR g (t)/dt\. Each panel in figure [2] shows 



time evolutions of T g (t) and Th(t) for < t < t npt . The black hole temperature at t npt is denoted as T g (t npt ) 
which is also attached in each panel in figure [2] The time t np t obtained from our numerical results are listed; 



Rh/Rg(0) 


87.775 


100 


10 3 


10 5 


10 6 


tnpt 1 tempty 


176.5 


2.039 


0.9989 


0.9155 


0.5363 


R g 


2.441 


2.549 


5.493 


25.49 


54.93 


Rh i tnpt 


1.2 x 10~ 7 


1.2 x 10~ 5 


2.4 x 10~ 4 


2.7 x 10~ 2 


4.6 x 10- 1 



(24) 



For the cases of Rh/R g (0) = 87.775, 100, 10 3 and 10 5 , the numerical plots stopped by v(t np t) = 1, but for the case 
of Rh/R g (0) — 10 6 , it stopped by Cx(t np t) + Cy^(i„ p4 ) = 0. This list (of t npt ) supports the discussion given at 
the beginning of this section UlI CI that the larger the radius Rh, the faster the black hole evaporation proceeds (the 
shorter the time t np t)- The third line in this list is for R* = l/(47rT*), which gives an important information in next 
subsection. The lowest line shows that our numerical results are consistent with the fast propagation assumption. To 
understand this statement, consider a typical time scale t ra d in which one particle of radiation fields travels in the 
hollow region from a black hole to a heat bath. The time scale t ra d can be given by t rac i — Rh, since the radiation 
fields are massless. Then, if the time t ra d is shorter than the time t npt , it is appropriate to consider that the fast 
propagation assumption is reasonable. Then the ratio t ra d/t npt (= Rh/tnpt) shown at the lowest line in list 
indicates that the fast propagation assumption is well satisfied. 
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FIG. 3: \C X + C^l - for the cases R h /R g (0) = 100, 10 3 and 10 5 in figure^] 
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FIG. 4: Jnpt/Jempty for each case in figure [2] The graph in upper-right panel is continuous and smooth. 



If the noncquilibrium region is ignored and the equilibrium model used in section [II] is considered with the same 
setting parameters, we find \C g \ = 2tt x 100 2 > 1000 = Ch, and the black hole settles down to a stable equilibrium 
state with a heat bath as explained in section [III That is, the occurrence of accelerated increase of temperature T g 
in figure [5] is obviously due to the nonequilibrium effect of energy exchange between black hole and heat bath. For 
the case Rh = 70i? g (0), the black hole settles down to a stable equilibrium state, since \C X +C ( Y 9) \ > C ( y l) holds 
for all time in this case (recall the discussion given in inequality (|14p ). As the radius Rh is set larger, inequality 
(fT4")) comes to be satisfied in the course of time evolution of the system. Then the black hole evaporation can be 
observed as an accelerated increase of temperature difference T g — Th for the case Rh = 87.775 R g (0). The time scale 
tnpt of this case is very longer than the time t empty- However, as the radius Rh is set larger and the nonequilibrium 
region becomes larger, the time scale t npt becomes shorter and we find t np t — t em pty about Rh — 10 3 R g (0). And t np t 
becomes shorter to t np t ~ 0.5363 t empty at Rh = 10 6 R g (0). If we set Rh > 10 7 , the combined heat capacity becomes 

positive Cx + Cy^ > which is forbidden in the framework of NPT model as explained in inequality (fT!?|) . That is, a 
black hole with nonequilibrium region of Rh > 10 7 can not be treated by the NPT model, and, as discussed in section 
IIIIBI (the paragraph one before the last one) , a black hole evaporation for this case should be described as a highly 
nonequilibrium dynamical process in which the black hole can not be treated with a stationary black hole solutions 
of the Einstein equation. The larger the nonequilibrium region, the faster the black hole evaporation process evolves 
into a highly nonequilibrium dynamical stage. 

Figure M shows time evolutions of \C X + C^ ] \ - f or the cases R h /R g (0) = 100, 10 3 and 10 5 . It is found from 
this figure that, as discussed at the end of section [III Bl inequality (| 14[) is not always satisfied during a black hole 
evaporation process. Even if inequality (|14[) is not satisfied at initial time, it must come to be satisfied in the course 
of evaporation process. When inequality (I14[) is not satisfied at initial time, the time t np t tends to be longer. 

Figure H] shows time evolutions of the ratio of energy emission rates, Jnpt/Jempty = (C g /Cx) (l — T^/T^), given 
by equation (|18p . It is recognised from this figure that the expectation given at the end of previous subsection is 
supported. Figure 2] indicates that, when a black hole evaporates, the relation J np t > Jempty comes to hold even if 
the converse relation J np t < Jempty holds at initial time. For the case Rh — 70i? g (0) in which a black hole stabilises 
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with a heat bath, the energy emission rate J npt disappears, J np t/ J empty — ► 0, as easily expected by the behaviour 
T g — * 7), =>■ J„ pf — * 0. 

5. Abrupt catastrophic evaporation 

In order to find a specific nonequilibrium phenomenon of NPT model, we point out a numerical evidence shown 
in list (|24|) that the black hole radius R* at time t npt is greater than unity R* > 1. According to the discussion in 
section UlI Bl (the paragraph one before the last one), when a black hole evaporates in the framework of NPT model, 
a highly nonequilibrium dynamical stage of evaporation process occurs at semi-classical level R* > 1. After the time 
tnpt, a black hole and radiation fields should be described as highly dynamical ones. 

In the following discussion, we make two steps: firstly, we show analytically that R* > 0(1) holds with a rather 
general condition. Secondly, a physical implication of R* > 0(1) is discussed. Then, a specific nonequilibrium 
phenomenon of NPT model is suggested. 

For the first, we analyse the energy transport equations ©■ Because of definition (|23p of time t npt , we consider 
two cases, t np t = t\ and t np t = t2, where t\ and t 2 are given in definition (|23p . But before proceeding to the analysis 
of these cases, we should point out the following: as explained at the end of appendix [Fj there is a possibility that 
equation Cx(R g ) + Cy (Rg) = has one, two or three solutions of R g for a sufficiently small Th- However, as 
discussed at the end of appendix [Fj even if there are two or three solutions of R g for our choice of Th, it is the lowest 
solution at which the quasi-equilibrium assumption is violated. Hence we consider the time t\ in definition (f2"3")) of 

t npt as the lowest solution of equation Cx(t) + Cy (t) = 0. 

Here we estimate the order of R* Consider the case t np t = h, where R* = R g (h) and Cx{R* g ) + Cy\R*) = 0. 
Because of Cy' > by definition, Cx{R* g ) — —Cy(R*) < holds. Consequently, according to a behaviour of 
Cx(Rg) explained in appendix [Cl (left panel in figure [7]), we find R* g > R g ~ 0.055 x (NR h ) 1/3 . Hence together with 
equation ([3]), we find R* > 0(1) for the situation Rh > O(10 3 ). And next consider the other case t np t = t 2 , where 

\Rg(t2)\ = 1. Because of t 2 <t%, we find Cx{t 2 ) + Cyfa) < 0. Then, because it is assumed that Th is small enough 
so that the validity condition (fT2| holds, we find by appendix [Fl (figure [8]) that R* = R g (t2) > R g o holds, where 

R g o is given by Cx(R g o) + Cy\R g o) = 0. Therefore, following the same discussion given for the case t np t — t\, we 
obtain R* > R g o > R g ~ 0(1) =>• R* > 0(1) for the situation Rh > 0(1O 3 ). In summary, the black hole radius R* 
at time t np t is greater than order unity R* > O(l), when the black hole evaporates in the framework of NPT model 
with the condition Rh > 0(1O 3 ). Here we have to note two points: first is that the condition Rh > 0(1O 3 ) is not a 
necessary condition but a sufficient condition for R* > 0(1), and there remains a possibility that R* > 0(1) holds 
even if Rh < 0(1O 3 ). Second point is that, if the nonequilibrium nature of black hole evaporation is not taken into 
account, the radius R* can not be greater than order unity but it becomes less than Planck length (see appendix [U| . 
The relation R* > 0(1) is peculiar to the NPT model. 

We proceed to the second part of this subsection, an implication of the above result R* > 0(1). Recall that 
a highly nonequilibrium dynamical stage of evaporation process occurs after the time t np t- (In the following, a 
"highly dynamical evaporation stage" means a "highly nonequilibrium dynamical stage of black hole evaporation 
process".) Because of R* > 0(1), a semi-classical (but not quasi-equilibrium) discussion can be available for the 
highly dynamical evaporation stage while the black hole radius shrinks from R* (> 0(1)) to R p (~ 1, Planck length). 
Then it is appropriate to consider that the mass energy of black hole evolves from E* (= R g /2) to E p (= R p /2 ~ 1, 
Planck energy), and that the speed of shrinkage of the horizon v — \dR g /dt\ is approximately unity v *~ 1 during the 
highly dynamical evaporation stage. The energy AE g = E* — E p is emitted during the highly dynamical evaporation 
stage. Further, for example, figure 2] of our numerical example in previous subsection implies a very strong luminosity 
Jnpt of Hawking radiation in the NPT model in comparison with the luminosity J e m P ty in the evaporation in an empty 
space. And note that the luminosity J em pty in an empty space is very strong as explained in section 1 of reference 
8]. Hence. J npt may be a huge luminosity. That is, the energy emission by a black hole in the NPT model may be 
understand as a strong "burst" . In addition to the luminosity of Hawking radiation, we consider the duration 5t of 
the highly dynamical evaporation stage. The duration St is estimated as 5t ~ R* g / V ~ Rg, and the following relation 
is obtained, 

5r= - ~ -«L < Ml = !_ « 0( o.oi) , (25) 

ternpty t empty t empty 12807ri? g (0) /! 400i? g (0)^ 

where we set N ~ 10 due to equation ([3]), and used R* < R g {0) and 1 <§; R g (0) (the initial radius R g (0) should 
be large enough to consider a semi-classical evaporation). This denotes St <C t empty . Further, for example, we find 
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St <C 0(0.01) < t np t/t em pty (=> St <C t npt ) for each panel in figure [2] Hence it seems reasonable to consider that St 
is very shorter than t npt . (For example it seems that the tangent dT g /dt seen in each panel in figure becomes very 
large quickly as t — > t npt , and T g will reach Planck temperature quickly just after t npt .) Hence it is concluded that 
the energy AE g = E* — E p bursts out of black hole with a very strong luminosity within St which is negligibly short 
in comparison with t npt . 

From the above, we can suggest the following picture: when a black hole evaporates in the framework of NPT 
model with a condition Rh > O(10 3 ), a quasi-equilibrium evaporation stage continues until t np t- Then a highly 
dynamical evaporation stage occurs at t npt . In the highly dynamical evaporation stage, a semi-classical black hole 
of radius R* > 0(1) evaporates abruptly (within a negligibly short time scale St) to become a quantum one. This 
abrupt evaporation in the highly dynamical evaporation stage is accompanied by a burst of energy A£" g . We call 
this phenomenon the abrupt catastrophic evaporation at semi-classical level R* > O(l), where "catastrophic" means 
that the speed of shrinkage of horizon is very high v ~ 1 and the energy AE g bursts out of black hole with a huge 
luminosity within a negligibly short time scale St in comparison with t npt . 

The above discussion is based on the NPT model. As shown in appendix[G] in the equilibrium model used in section 
ITTl and in the model that a black hole evaporates in an empty space, the black hole radius becomes Planck size before 
the shrinkage speed of horizon reaches unity. The highly dynamical evaporation stage and the abrupt catastrophic 
evaporation at semi-classical level do not occur in both models. Hence the abrupt catastrophic evaporation at semi- 
classical level seems to be a specific nonequilibrium phenomenon of NPT model. 

Finally we estimate a typical time scale of black hole evaporation in a heat bath. It is reasonable to consider 
the duration of quantum evaporation stage is about Planck time. Then, a time scale of evaporation is given by 
t n pt + St + O(l) ~ t npt . The time t npt gives a typical time scale of black hole evaporation in a heat bath. 

6. Conclusion of section WTC] 

From the above, we find the nonequilibrium effect of energy exchange on a black hole evaporation in a heat bath 
as follows: When the nonequilibrium region around black hole is not so large, the evaporation time scale t npt in a 
heat bath becomes longer than the evaporation time scale t emp t y in an empty space ( ignoring gravitational interactions 
between black hole and radiation fields). However, as the nonequilibrium region around black hole is set larger, the 
time scale t np t becomes shorter than the time scale t emp t y ■ That is, the nonequilibrium effect of energy exchange tends 
to accelerate the black hole evaporation. Further, for the case Rh > 0(1O 3 ), a quasi-equilibrium evaporation process 
proceeds to a quantum evaporation stage abruptly at t npt at semi- classical level R* > 0(1) with a burst of energy (the 
abrupt catastrophic evaporation). 

Finally in this section, we discuss about a black hole evaporation in an empty space in a full general relativistic 
framework. Note that, even if a black hole is in an empty space, there should exists an incoming energy flow into 
the black hole due to the curvature scattering. When the curvature scattering is taken into account for the case of a 
black hole evaporation in an empty space, we can interpret the whole system as that a black hole is surrounded by 
some nonequilibrium matter fields which possess outgoing and incoming energy flows of Hawking radiation under the 
effects of curvature scattering. Further, since the curvature scattering occurs whole over the spacetime, it is expected 
that the nonequilibrium region is so large that a condition which corresponds to Rh > O(10 3 ) in the NPT model 
holds. Hence, if the NPT model is extended to a full general relativistic model, we can expect that a black hole 
evaporation in an empty space can be treated in the framework of full general relativistic version of the NPT model 
(with removing the heat bath), and that the abrupt catastrophic evaporation at semi-classical level may occur as well 
since the nonequilibrium region is sufficiently large. 

D. Final fate of black hole evaporation 

In previous section, the black hole evaporation processes before a quantum evaporation stage have been considered. 
In this section, we consider whether some remnant remains after the quantum evaporation stage or not. The study 
on the final fate of black hole evaporation, which has a close relation with the information loss problem [j| , is usually 
carried out in the context of quantum gravity. However, since no complete theory of quantum gravity has yet been 
constructed, it is meaningful to some extent to study the final fate of black hole evaporation with an appropriate 
model of black hole evaporation without referring to present incomplete theories of quantum gravity. Therefore, we 
use the NPT model and try to extract what we can suggest about the final fate of black hole evaporation. 

For the time being, we assume that a black hole evaporates out completely and only equilibrium radiation fields 
remain as the end state of quantum evaporation stage. If a contradiction results from this assumption, we may 
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FIG. 5: If a black hole evaporates out completely, when a black hole temperature T g reaches Tg, the black hole is suddenly 
replaced by radiation fields of temperature T* ad which is determined by the energy conservation. However it is impossible from 
entropic viewpoint. 



conclude that a remnant will remain after a black hole evaporation. This section aims to suggest a necessity of some 
remnant by the reductive absurdity. 

As discussed at the end of section ITlI CI it can be expected that the abrupt catastrophic evaporation at semi-classical 
level occurs in a full general relativistic framework. Therefore we consider the case Rh > O(10 3 ) in the NPT model 
which denotes the occurrence of abrupt catastrophic evaporation at semi-classical level. Then, under the assumption 
of complete evaporation of black hole, we can draw a scenario in the framework of NPT model as follows. 

As a quasi-equilibrium stage of black hole evaporation proceeds until t npt , the temperature T g approaches a critical 
value T* = l/47ri?*, where R* = R g (t np t). Then, a highly dynamical evaporation stage and a quantum evaporation 
stage follow. Under the assumption of complete evaporation, these successive stages are described as one phenomenon; 
a black hole emits completely its mass energy E* (= R g /2) within a negligibly short time scale just after the time 
tnpt (abrupt catastrophic evaporation), and equilibrium radiation fields remain. In the NPT model, this abrupt 
catastrophic evaporation under the assumption of complete evaporation is treated by the procedure; replace a black 
hole of radius R* with equilibrium radiation fields of volume V* = (47r/3) R g 3 - In this replacement, thermodynamic 
states of heat bath and nonequilibrium radiation fields around the region of V* are not changed (see left panel in 
figure [5]) . Here note that the energy of equilibrium radiation fields in volume V* should equal the mass energy E* of 
black hole. Therefore the temperature T* ad of equilibrium radiation fields is determined by 

/ 3 \ 1/4 

E*=4a'T^ d V* => T; ad = — — — , , (26) 



32a'irR 



g ' 

where 4c' T 4 V is the equilibrium energy of radiation fields (see appendix [B]) . When a sufficiently long time has 
passed after the abrupt catastrophic evaporation (under the assumption of complete evaporation), the whole system 
reaches an equilibrium state in which a heat bath and radiation fields in the hollow region have the same equilibrium 
temperature. However, without considering this totally equilibrium end state of the whole system, but with considering 
the states just before and just after the abrupt catastrophic evaporation, we can discuss whether a remnant remains 
or not after a black hole evaporation. 

Here we should notice that the total entropy must increase along a black hole evaporation process, since the whole 
system which consists of a black hole, a heat bath and radiation fields is isolated from outside environments. Therefore 
the relation S£ ot < S*^ should hold, where «S t * ot is the total entropy of the whole system at time t np t just before the 
abrupt catastrophic evaporation, and SfJ t is the total entropy of the whole system just after the complete evaporation 
of black hole. Because the abrupt catastrophic evaporation is described by a simple replacement of a black hole with 
equilibrium radiation fields, we find that the entropy difference AS% ot = S^ t — S^ ot is given by 

^^tot = '-'tot ~ "tot — '-'rad — D g I \ z 1 ) 

where S* is the black hole entropy of temperature T* and S* ad is the equilibrium entropy of radiation fields of volume 
V* and temperature T* ad . If the complete evaporation of black hole is true of the case, a relation AS^ ot > must 
hold. However it is impossible as shown below. 

Using equations of states ([1]) and the equilibrium entropy of radiation fields (16cr'/3) T 3 V together with the tem- 
perature T* ad of equation (f2"6")l , the entropy difference is calculated as 

As; ot = T ; a \ v; -s; = \( ^ ) ^ r; 3 / 2 - * r; 2 . (28) 
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A schematic graph of AS^ ot is shown in right panel in figure [5j where R* = (8/27-\/57r) y/N ~ 0.075 x y/~N is given 
by AS; ot = and a' = Nn 2 /120. We find AS? ot > for R* < R* and AS; ot < for R* > R* g0 . Here recall that, 

because of Cy^ > by definition and a behaviour of Cx(R g ) explained in appendix [C] (left panel in figure [7|), an 
inequality R* > R g holds due to the validity condition (jl2[) . where R g is given by solving Cx = for R g . Further, 
using equation (jUTTI) . we find R g /R* ~ (27^/5^/32(30^) x / 3 ) x (R 2 /N)^ 6 ~ 0.734 x {R 2 /N)^ 6 . Then, by equation 
(|3"|) and requirement Rh > O(10 3 ) of the occurrence of abrupt catastrophic evaporation, we obtain 

§->W= Oil)- (29) 

This indicates AS* ot < 0. Consequently, by the reductive absurdity as explained at the beginning of this section, a 
complete evaporation of black hole is impossible. 

From the above, we give a suggestion as follows: a complete evaporation of black hole is prohibited and a remnant 
should remain after a black hole evaporation. The entropy of remnant should guarantee the increase of total entropy. 

We comment about the mass of remnant. One may think that the mass energy (or internal energy) of remnant may 
be extracted by applying an energetic analysis to the "entropic" discussion given above. One may obtain a minimum 
entropy S^™^ of remnant as <S^„ = S g , and a mass energy of remnant may be obtained from S^™'. However, 
because equations of states of remnant is unknown, it is impossible to obtain an energy of remnant from the entropic 
discussion. On the other hand, as discussed in section [III CI (before equation [)25p ). the mass of quantum black hole is 
expected to be about one Planck energy. Therefore, as long as the NPT model is extrapolated over its validity, the 
mass of remnant seems to be of the order of one Planck energy. 

In the rest of this section, we discuss more about the validity of discussion given in deriving inequality (|29p . There 
are four points we are going to discuss here. For the first, we discuss about the entropic viewpoint in NPT model. As 
explained in appendix the NPT model is constructed referring to the validity of energetics of equilibrium model 
used in section [H] and reference [l| . Therefore one may think that the NPT model can not reflect the entropy of 
gravitational field around a black hole. However, as mentioned in equation ([2|), using the Euclidean path integral 
method for a black hole spacetime [ic| , the entropy of gravitational field itself is just given by the quantity S g appeared 
in equations iJTJ. That is, the black hole entropy is the entropy of the whole gravitational field in a black hole space 
time. Hence it seems that an entropic discussion in NPT model is appropriate. 

Next for the second point, we discuss about the usage of NPT model in this section. If we want to analyse details 
of the final stage of black hole evaporation process, then a quantum gravity is necessary, since a black hole becomes 
quantum, R g < 1. On the other hand, because of the quasi-equilibrium assumption, the black hole radius has to be 
restricted as (R g (0) >)R g > 0(1) in the framework of NPT model. However this restriction R* > 0(1) does not 
mean that the NPT model is never available for considering the final fate of black hole evaporation. Although the 
NPT model is based on a classical Schwarzschild black hole, once one assume that a black hole will evaporate out 
completely, then it is true under this assumption that radiation fields remain and its entropy can be counted after a 
quantum evaporation process ended. That is, even though we never refer to any details of present incomplete theories 
of quantum gravity, we can compare an entropy of black hole before the start of quantum evaporation process with 
an entropy of radiation fields after the end of quantum evaporation process. Hence our analysis done above seems to 
be reasonable to approach the final fate of black hole evaporation. 

For the third point, we turn our discussion to the starting assumption that a black hole evaporates out completely 
and equilibrium radiation fields remain. As seen above, equation (|28p has led a contradiction to deny the assumption 
and to result in necessity of a remnant. Here recall that equation (|28| is calculated with using the equilibrium entropy 
of radiation fields. Then, one may think that it is more general to modify the starting assumption so that radiation 
fields after the complete evaporation of black hole are not necessarily in an equilibrium state. However, as explained 
in detail in reference [l2j], the nonequilibrium entropy of radiation fields should be smaller than the equilibrium one 
S* ad (< S*) used above, since the equilibrium state is the maximum entropy state. That is, it is reasonable to 
expect that an inequality AS^ ot < holds stronger for a modified starting assumption which requires nonequilibrium 
radiation fields after a complete evaporation of black hole. 

Finally we discuss about a black hole evaporation in an empty space (a situation without heat bath). As discussed 
at the end of section Ull C[ the abrupt catastrophic evaporation at semi-classical level is expected to occur for a black 
hole evaporation in an empty space in a full general relativistic framework. Further note that calculations of entropy 
difference AS^ ot in equation (|28p do not depend on the outside of black hole, but depend only on the black hole entropy 
S* and the equilibrium entropy S* ad of radiation fields. That is, without respect to the degree of nonequilibrium 
nature of the environment around black hole, once the abrupt catastrophic evaporation of black hole occurs, the same 
discussion as given in this section seems to be applicable to a black hole evaporation in an empty space. Therefore, 
a full general relativistic treatment of black hole evaporation in an empty space may result in necessity of a remnant 
at the end state of black hole evaporation. 
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IV. SUMMARY AND DISCUSSIONS 

When a black hole evaporates in a heat bath, there exists an energy flow between the black hole and the heat bath 
due to the Hawking radiation and the accretion of matters from heat bath to black hole. This denotes that the whole 
system which consists of a black hole and a heat bath is not in an equilibrium. We have investigated nonequilibrium 
effects of energy flow on the black hole phase transition. Our analysis is carried out with the NPT model which 
reflects the nonequilibrium nature of energy exchange between a black hole and a heat bath. Then it is found that, 
contrary to the original equilibrium model used in reference [l[ , inequality (|14p evaluated at initial time can not work 
as a criterion of the occurrence of instability of equilibrium state of the whole system. The nonequilibrium effect of 
energy exchange tends to make the situation more dynamical and complicated. 

For the case that a black hole evaporates in a heat bath, the energetics of NPT model gives an understanding of 
nonequilibrium nature of black hole evaporation as follows: 

• When the nonequilibrium region around black hole is not so large, the evaporation time scale t np t in a heat bath 
becomes longer than the evaporation time scale t empt y in an empty space, because of the incoming energy flow 
from the heat bath to the black hole. However, as the nonequilibrium region around black hole is set larger, 
the time scale t np t becomes shorter than the time scale t emp t y , because a nonequilibrium effect of temperature 
difference between the black hole and the heat bath appears as a strong energy extraction from the black hole 
by the heat bath. That is, the nonequilibrium effect of energy exchange between the black hole and the heat 
bath tends to accelerate the black hole evaporation. Further for the case Rh > O(10 3 ), the abrupt catastrophic 
evaporation at semi-classical level occurs. 

Concerning a black hole evaporation in an empty space in a full general relativistic framework, we can expect this 
understanding holds due to the curvature scattering, as discussed at the end of section |HI CI 

Two possible future applications of our discussion are suggested here. First possible future application is primordial 
black holes formed in radiation fields in the early universe [3| • The lower bound of the mass of primordial black holes 
which survive at present is estimated usually using the time scale t emp t y of black hole evaporation in an empty space. 
However, since the evaporation time scale t np t in a heat bath changes according to how the radiation field in the early 
universe is far from an equilibrium, the lower bound of mass of primordial black holes may not be determined simply. 
Second possible future application is black holes created in particle accelerators in the context of TeV gravity j4j. 
The black holes in particle accelerators, if they exist, will be created in a quark-gluon plasma in our brane and in a 
graviton gas in the higher dimensional spacetime. It is usually argued in many papers that, although a black hole in 
accelerator is highly dynamical at the moment of its creation, it will settle down to a stationary black hole before it 
proceeds to a quantum evaporation stage, since the evaporation time scale t emp t v in an empty space is slightly longer 
than the quasi-normal ringing time scale. However, since the evaporation time scale t npt in a heat bath changes 
according to how the quark-gluon plasma and the graviton gas are far from an equilibrium at the moment of particle 
collision, the black holes in particle accelerators may not settle down to a stationary black hole before its quantum 
evaporation stage. 

If a black hole is in an expanding universe as primordial black holes, not only the energy exchange between a 
black hole and its environments but also the cosmological expansion may cause nonequilibrium effects on black hole 
evaporation. Then one may think that the evaporation process may become a highly nonequilibrium one and be 
accelerated faster than the case without cosmological expansion. It is true that cosmological expansion disturbs the 
thermal nature of Hawking radiation and the spectrum of energy flow from black hole can not be treated as a Planckian 
distribution [l3| . However nonequilibrium effects of cosmological expansion should require different considerations 
from the idea of this paper. Time evolution of the acceleration of cosmological expansion may be essential to determine 
whether the energy emission rate by black hole is enhanced or weakened. If the time derivative of acceleration of 
cosmological expansion is positive (or negative), the emission rate may be weakened (or enhanced) That is, 

when one consider details of black hole evaporation, the origin of nonequilibrium nature may be important. It seems 
that there are a variety of nonequilibrium effects. It is the nonequilibrium effect of energy exchange (or temperature 
difference) between a black hole and its environments that we investigates in this paper. 

Next, concerning the end state of black hole evaporation, the entropic discussion of NPT model gives the suggestion 
as follow: 

• Without respect to the size of nonequilibrium region around black hole, a remnant should remain after a black 
hole evaporation. The entropy of remnant should guarantee the increase of total entropy. 

As long as the NPT model is extrapolated to a highly nonequilibrium dynamical stage and a quantum stage of 
evaporation after the time t np t, the mass of remnant seems to be of the order of one Planck energy. And, concerning 
a black hole evaporation in an empty space in a full general relativistic framework, we can expect the suggestion of 
necessity of remnant holds due to the curvature scattering, as discussed at the end of section UlI Dl Further, recalling 
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the third paragraph of section [IJ our suggestion implies that a remnant can preserve the initial information of black 
hole formation, and that the unitary invariance of quantum theory is not violated. That is, the information loss 
problem [f| may not exist. 

Finally we have to comment on another possibility ("non-remnant" possibility) to preserve the unitary invariance 
along a black hole evaporation. The detailed properties of Hawking radiation may preserve the initial condition of 
gravitational collapse. Because it is the observation at an infinite future that a thermal radiation is observed in 
the Hawking radiation Q, a deviation of observed spectrum from Planckian distribution (at finite future from the 
onset of gravitational collapse) may carry the information of initial condition into somewhere on the spacetime, and 
the unitarity may be retained over the whole spacetime. However, if a naive and simple extrapolation of the NPT 
model to a quantum evaporation stage is valid (as assumed in section fill Dp . then the thermodynamic state of outside 
environment of black hole is composed of a steady state of radiation fields and an equilibrium state of heat bath. 
Hence, because the steady state is described by a simple superposition of equilibrium states of different temperatures 
(see appendix [B]) , it seems that a complicated (highly noncquilibrium) information of initial condition can not be 
preserved in steady states of radiation fields and equilibrium states of heat bath. Therefore, as far as the NPT 
model is extrapolated over its validity, a possible conclusion of A5 t * ot < (see equation (|29|) ) is not a complicated 
nonequilibrium state of outside environment but a remnant. 

APPENDIX A: ON GRAVITATIONAL INTERACTIONS BETWEEN A BLACK HOLE AND A MATTER 

FIELD 

This appendix tries to give an intuitive explanation why a primitive model considered in section |TT] and reference 
[l| retains the essence of black hole phase transition. Since section [TT1 and reference [l[ do not include any relativistic 
effects except for equations of states ([1]) , it may be expected that the essence of black hole phase transition is roughly 
explained with a Newtonian discussion. 

But before proceeding to a Newtonian discussion, we consider the Hawking radiation by a Schwarzschild black hole 
in the context of general relativity. When a monochromatic wave mode of matter field of Hawking radiation propagates 
from a point near the black hole horizon to asymptotically flat region, a gravitational Doppler effect causes a red shift 
on the mode, uj t = cj/^/l — r g /r, where lo t and u> are respectively frequencies at starting point and at asymptotically 
flat region, r g is the gravitational radius of black hole, and r is the radial coordinate (areal radius) at starting point. 
It is obvious lo t — ► oo as r — ► r g . It has already been known that [ToT ]. when a black hole is in a heat bath and they 
are in an equilibrium, the local equilibrium temperature T(r) of heat bath at a point of radial coordinate r should be 
given by T(r) = T g / 'wl — r g /r, where T g is a usual Hawking temperature of black hole which appears in equations 
(Q]), and the factor yl — r g /r is due to the gravitational red shift. We can understand T(r) as that, in order to retain 
a thermal equilibrium against attractive gravitational force, a temperature which is higher than the asymptotic value 
should be required, since the higher temperature denotes the higher pressure against gravitational force. 

Next we shift to a Newtonian discussion about matter field which composes a heat bath around black hole. Consider 
the case that the matter field is in an equilibrium with black hole. We put three assumptions: first is that the whole 
system which consists of the black hole and the matter field is isolated (total energy is conserved as the primitive 
model used in [l[). Second assumption is that, for simplicity, the matter field of heat bath is the same with that 
of Hawking radiation. Third one is that self-interaction among particles of the matter field is negligibly weak in 
comparison with particle's kinetic energy and gravitational interaction with black hole. 

Because of the assumption of very weak self-interaction, the energy e of one particle of the matter field is 

e = K(r) + U g (r), (Al) 

where K(r) is kinetic energy at radial distance r from the center of black hole, and U g (r) is gravitational potential 
given by black hole. According to the gravitational Doppler effect tu r = lu / ^Jl — r g /r, the gravitational potential 
should have the properties, U g — > — oo as r — ► r g , and U g — > as r — > oo. Then the total energy of the whole system 
is given by E to t = E g + Ee, where E g is black hole mass, and a summation £ is taken for all composite particles of 
matter field. Here note that, since the whole system is in an equilibrium, a black hole is static and its mass is constant, 
Eg = constatnt. Further, because of the assumption that the whole system is isolated, we have E tot = constant. 
Therefore we find £e = constant, and together with the assumption of very weak interaction (energy is not exchanged 
among the particles), it is implied e = constant. Then, due to the asymptotic value of U g (r — > oo) = 0, we find 
e — K(r) + U g (r) — Koa, where = K(r — > oo) is the value of kinetic energy at infinity. 

In order to understand the meaning of ifoo, we consider the local temperature of matter field around black hole. 
Due to the law of equipartition of energy, we have 



(K(r)) ocr(r), 



(A2) 
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where (K(r)) is a statistical mean value of K(r) at r, and T(r) is the local equilibrium temperature of matter field 
at r. Then, because of relation K(r) = e — U g (r) = — U g (r), we find properties of local equilibrium temperature, 
T(r) — > oo as r — > r g and T(r) — > T g as r — > oo. These behaviours of T(r) correspond to those of a general 
relativistic relation T(r) = T g /wl — r g /r. This may be interpreted as follows: when we consider separately a density 
of internal energy e(r) of matter field (where e(r) — n(r) (K(r)), n(r) is the number density of particle at r) and a 
density of gravitational interaction energy u(r) by black hole (where u(r) — n(r) U g (r)), then infinite values appear 
e(r — » r g ) = oo and u(r — ► r g ) = — oo, and one may think thermodynamic considerations break down at r = r g 
due to the infinite values. However, the relation K(r) + U g {r) = Koo implies that these infinite values cancel each 
other. That is, although a realistic situation of matter field around black hole should be understood that its energy 
is influenced by the gravitational potential by black hole; however, when one considers an "effective" internal energy 
of matter field Eoo = N (Koo) (where N is the number of particles of the matter field), then this energy Eoo equals 
the sum of the whole internal energy of the matter field and the interaction between black hole and particles of the 
matter field, Se = Eoo- That is, the infinite values of energies at r — r g are "renormalised" into the effective value 
Eoo. This "renormalised energy" Eoo is the internal energy of the matter field evaluated by a thermal equilibrium 
state of asymptotic temperature T g without the influence of gravitational field of black hole. Hence, the total energy 
is calculated as E tot = E g + E e = E g + E^ . 

From the above it is recognised that, as long as we use the renormalised internal energy Eoo of matter field of 
asymptotic temperature T g , energetics of the whole system which consists of a black hole and matter field can be 
considered with ignoring the gravitational interaction between black hole and matter field. Here recall that section 
HT1 and reference |l| consider the energetics of equilibrium model in which a black body of temperature T g is put in a 
heat bath of the same temperature T g and gravitational interaction between them is ignored. Hence section [TTI and 
reference [l[ can retain the essence of black hole phase transition. It is the self-gravitational effect of black hole on its 
own thermodynamic state that causes the black hole phase transition. 

So far we have considered an equilibrium between a black hole and a matter field. Finally turn our discussion 
to the case that the matter field is not in an equilibrium with the black hole. In the nonequilibrium case it is not 
generally clarified whether equation (|A2|) holds or not, since nonequilibrium statistical mechanics has not yet been 
established in a general form. However when we consider non-self-interacting massless matter fields (named as the 
radiation fields in the NPT model), then we can treat the nonequilibrium effect of matter fields by two-temperature 
steady state thermodynamics for the radiation fields which is summarised in appendix [B] Hence it is expected 
that the original model considered in reference [l| can be extended to include the nonequilibrium effects of Hawking 
radiation with utilising the steady state thermodynamics for the radiation fields. 

APPENDIX B: STEADY STATE THERMODYNAMICS FOR THE RADIATION FIELDS 

This appendix summarises two-temperature steady state thermodynamics for a radiation field fl2| . As done in refer- 
ence [lH, we consider a photon gas as the radiation field. However it is easily extended to general non-self-interacting 
massless matter fields by replacing the Stefan-Boltzmann constant a with the generalised one a' as mentioned at the 
end of section Q] 

We consider the following model which is named SST after the Steady State Thermodynamics. 

SST model: Make a vacuum cavity in a large black body of equilibrium temperature T out and put an another smaller 
black body of temperature Ti n T out ) in the cavity (see figure [5]). For the case Tj„ > T out , the radiation field 
(photon gas) emitted by two black bodies causes an energy flow from the inner black body to the outer one. 
When the outer black body is isolated from outside environments, the whole system which consists of two black 
bodies and radiation field between them relaxes to a total equilibrium state in which two black bodies and 
radiation field have the same equilibrium temperature. 

Note that equations of states of two black bodies are not specified, and we assume that gravitational interactions 
among two bodies and radiation field is negligible. Further we consider the case satisfying the following two conditions: 
the first is that thermodynamic state of each black body changes along a sequence of equilibrium states during the 
relaxation process of the whole system, and the second condition is that the volume of cavity is so small that the speed 
of light is approximated as infinity during the relaxation process. Then, as shown in reference thermodynamic 
state of radiation field sandwiched by two black bodies should change along a sequence of nonequilibrium states 
during the relaxation process. Further due to the Stefan-Boltzmann law, each nonequilibrium state in the sequence 
of radiation field's states has an energy flow J sst — a ( T£ — T^ ut ) A in , where A in is the surface area of inner black 
body. This J ss t equals the net energy exchanged par a unit time between the two black bodies via the radiation field. 

From the above discussions, it is enough for us to consider a set (state space) of steady states (macroscopically 
stationary nonequilibrium state) of the radiation field whose nonequilibrium nature is raised by the steady energy 
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FIG. 6: SST model. Radiation field between outer and inner black bodies is in a steady state. 

flow J sst . That is, the sequence of nonequilibrium states of the radiation field during the relaxation process lies in 
the set (state space) of steady states. A consistent thermodynamic framework for the steady states of radiation field 
has already been constructed (T3 |. The outline of the construction of internal energy and entropy are as follows. 
For the first we consider a general bosonic gas. The energy Eb of the gas is given by, 

/dp 3 
(2 7r) 3 ^ X e ^' S ' 

where p is a momentum of particle, x is a spatial point, ep t s is an energy of particle of p at x, and gp^ and fp t s 
are respectively the number of states and the average number of particles at a point (p, x) in the phase space of the 
bosonic gas. Further, refer to the Landau-Lifshitz type definition of nonequilibrium entropy Sb for bosonic gas [lq . 



S h = 



f dp 3 



Note that it has also been shown in reference [Ty| that the maximisation of Sf, for an isolated system (SSf, = 0) 
gives the equilibrium Bose distribution. This is frequently referred in many works on nonequilibrium systems as the 
H-theorem. However in reference [l6j], concrete forms of gp y g and fp t $ are not specified, since an arbitrary system is 
considered. 

In the SST model, the radiation field is sandwiched by two black bodies, and the particle (photon) is massless and 
collisionless. Then we can determine ep t s, gp,s and fp t g as 

tp.s = w , gp. s = 2 (helicities of a photon gas) , fp S = — — , (B3) 

exp[w/T{p,x)]-l 

where a frequency of photon u> = \p\, and T(p, x) is given by 

rp/^ * \Tin for p = p in at x 
[T out tor p = p out at a; 

where pi„ is a momentum of photon emitted by the inner black body and p Q ut by the outer black body. The in- 
dependence of T(p, x) arises, because the directions in which photons of pi n and p ou t can come to a point x vary from 
point to point. 

From the above, we obtain the steady state energy Eb and entropy Sb, 

E b = J dx 3 e rad (x) , e b (x) = Aa (g in (x) T^ n + g out {x) T* ut ) , (B5) 

and 

Sb = J dx 3 s rad {x) , s b (x) = (g in {x) T 3 n + g out {x) T 3 ut ) , (B6) 

where gi n (x) is a solid angle (divided by 47r) covered by directions of p lin at x as shown in figure [6l and g ut{x) is 
similarly defined with p ou t- By definition, we have gi n (x) + g ou t(x) = 1< Here it should be noted that, if the radiation 
field is in an equilibrium of temperature T and volume V , the energy E — AaT V and the entropy S = (16<t/3) T 3 V. 
Hence we find that, because photons are collisionless (non-self-interacting matter field), the steady state energy and 
entropy are given by a simple linear combination of the values which are calculated if the radiation field is in an 
equilibrium of temperature Ti„ or T out . 

Further, with defining the other variables like free energy with somewhat careful discussions, it has already been 
checked that the 0th, 1st, 2nd and 3rd laws of ordinary equilibrium thermodynamics are extended to include the steady 



20 



states of a radiation field. That is, the steady state thermodynamics for a radiation field has already been constructed 
in a consistent way. Especially on the steady state entropy (|B6|) . it has already been revealed in [12| that the total 
entropy of the whole system which consists of two black bodies and radiation filed increases monotonously during the 
relaxation process of the whole system, dSi n + dS ou t + dSb > 0, where the equality holds for total equilibrium states, 
Si n and S ou t are entropies of two black bodies which are given by ordinary equilibrium thermodynamics, and it is 
assumed that the two black bodies are made of ordinary materials and have positive heat capacities. 

When we consider a massless fermionic field instead of the radiation field in SST model, the formula of energy (|Bip 
holds for fermions as well. But the formula of entropy (|B2|) is replaced by [l6| 



f dp 3 

Sf = - J ( 27 r)3 dx3 9,y ' s t ^ 1 ~~ f p ' s ^ ln ^ 1 ~~ ^ + ln fe< s ^ ' ( ' B7 ' ) 

Then, following the same procedure as for a radiation field with using the distribution function of massless fermions 
fp,s — [ e ' XJ p[ UJ /T(p, x)] + 1 ] , we can obtain the steady state energy and entropy of a massless fermionic field as 

Ef = 8 ~~2 Eb ' Sf = 8 ~2 Sb ' < ^ B8 ' ) 

where gp t g — rif is the effective number of hclicities of the fermions. Hence, when the black bodies in SST model emit 
rif massless fermionic modes and n& massless bosonic modes (rib — 2 for a radiation field), the steady state energy 
and entropy of matter fields between two black bodies is given by equations (|B5[) and (|B6|) with replacing a by a' 
given at the end of section |TJ 

Finally for the convenience to follow equations (|B5|) . (|B6[) and (|B8|) . we list useful formulae, 

36 ' 

da; = , / da; \n(e x + l) = \- F , / dx lnfl + e x ) = F, 

eHl 120' Jo e x + l y ' 60 ' 7 1 + e~ x y ' 180 

where F = [(ln2) 2 - tt 2 ] (ln2) 2 /6 + 4 0(4, 1/2) + (71n2/2) C(3), ((z) is the zeta function and cj)(z, s) is the modified 
zeta function (Apell's function). 

APPENDIX C: HEAT CAPACITY C X (R S ) AS A FUNCTION OF R g 

This appendix shows a behaviour of heat capacity Cx of sub-system X as a function of black hole radius R g . The 
Cx(Rg) is given by equations @ and (fTT|) as 

C x (R g )=C g {R g ) + C^ d (R g ), (CI) 

where < R g < Rh by definition and 



' , x 3 n A f°° , x 2 , . B Utt 4 /- 00 a; 2 . ^ _„ 

da; = — , / da; ln(e — 1) = , / dx lnfl — e ) = 

e x - 1 15 ' .L e x -I K ' 180 ' ./„ 1 - e-* 1 ; 



C g (R g ) = -2wR 2 q 



J g\ ±l -g) 
<g) 

J rad 



(C2) 



By equations ([T]) and (O, this is rewritten into the following form, 



where 



Then we find 



C x (x)=-27rRlx 2 + ^f(x) , *=|jj, (C3) 



m = ±± \l-x 3 -(l-x 2 ) 3/2 ]+- (x~^ VA . (C4) 
3r L J a; V / 



Cx — > oo as x — > 



a' 



C x = -ZkRI + — — — at x = 1 



(C5) 
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c x (Rg) 



C x (R h ) = C,^(R h ) + constant 



R g * 0.055 (NR h )' 
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R„ 




FIG. 7: Left panel is for Cx(R g ) as a function of _R g , where 2s/2/3 ~ 0.943. Right is for Cx(Rh) as a function of where 
3^2/4 ~ 1.06. Due to the validity of NPT model, a' = O(l) and O(l) < i? 9 < R h , then Cx(-R 9 = < holds. 



Here recall the facts; Rh > R g by definition, (Rh >)R g > 0(1) due to the quasi-equilibrium assumption, and 
N = O(10) by equation ([3|). Then we can expect that an inequality Cx(x = 1) < holds for the NPT model. Further 
the differential of Cx (%) is given by 



dC x (x) 
dx 



-4ttRI x 



a' df(x) 



where 



df(x) 4-3x 2 ~4y/l-x 2 



8tt 2 dx ' 



df (x) „ 2^2 

• M ' < for x < 



dx 



x 4 VT 



dx 
df(x) 

dx 



> for x > 



3^ 
2V2 ' 



(C6) 



(C7) 



where 2^/2/3 ~ 0.943. Adjusting the values of Rh and N in a', an inequality dCx/dx > may hold for x > 2y/2/3. 
However, because of Cx(x = 1) < 0, it is obvious that the equation Cx = has only one solution. Hereafter 
R g (= Rh x) denotes the solution of Cx(R g ) = 0. A schematic graph of Cx(R g ) is shown in figure [7] (left panel). 

Finally in this appendix, we estimate the value of x. Since x < 1 by definition, we apply the Tayler expansion, 
(I — x 2 ) a = I — a x 2 + 0(x A ), to equation (|C3[) and obtain 



C x (x) = -2ttR 2 



(C8) 



where e = a'/l6iT 3 R 2 = N/1920ttR 2 . By equation © and R h > R g > 0(1) due to the quasi-equilibrium assumption, 
we find e < 1. Then Cx(x) becomes 



1 



-2wRi- (x 6 -e + eO(x) 
x 



C x (x) 

Hence, by taking leading terms of x and e, we find an approximate expression for x as 

C x (x) =0 

This gives an approximate value of R g , 



a/3 



Ra 



NR h 
1920tt 



1/3 



0.055 x (NR h 



,V3 



(C9) 



(CIO) 



(Cll) 



APPENDIX D: HEAT CAPACITY Cx(Rh) AS A FUNCTION OF R h 

This appendix shows a behaviour of heat capacity Cx of sub-system X as a function of outermost radius Rh of 
the hollow region in NPT model. We make use of the calculations done in appendix [Cl According to equation (|C3[) , 
Cx(Rh) is expressed as 



C x (R g ) = -2nRl + Cl 3 J d (R h ) 



8tt 



a fa) ■ 



(Dl) 
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where x = R g /Rh, < x < 1 by definition, and f(x) is given by equation (|C4[) . This denotes that Cx(Rh) behaves 
as C^J^Rh) + constant. The differential becomes 



Using equation (|C7|) , we find 



dCx(Rh) 
dR h 



^_ x 2 d f(x) 

8n dx 



(D2) 



(dC x _ a' 
~dR^ ~ 8^ 
dC x 

- dRh 



at Rh = R g [x = 1) 
oo as Rh — ► oo (x — ► 0) 



(D3) 



Hence referring to the limit values (|C5j) and the behaviour (|C7[) of df/dx, a schematic graph of Cx(Rh) is obtained 
as shown in figure [7] (right panel). This denotes that, since Cx < must hold by the validity of NPT model (fT2"| , the 
absolute value \Cx(Rh)\ is monotone increasing for Rh > (3v / 2/4)_R g 



1.06i? g in the framework of NPT model. 



APPENDIX E: PROOF OF THE INEQUALITY C g /Cx > 1 

This appendix gives a proof of inequality C g /Cx > 1 under the condition Cx < (see the discussion on inequality 
(fT2")l ). We make use of the calculations done in appendix [Cl By equation (|C3|) and definition of Cx given in equations 
(fTTj) , the heat capacity C!fJ d of out-going radiation fields is expressed as 

= -^m, (ei) 

where x = R g /Rh, < x < 1 by definition, and /(x) is given by equation (|C4jl . Then equation (|C7j) indicates 

C^(x) > C^(2V2/3). Here equation jC4} gives /(2a/2/3) ~ 0.845 =>■ /(2V2/3) > 0. Therefore we find > 
for < x < 1, which is consistent with a naive expectation that an ordinary matter like radiation fields has generally 
a positive heat capacity. Consequently the condition Cx < gives an inequality C g /Cx > 1 for < x < 1 as follows 
(see equation (J2]) for C g ), 

C x =C g +C { r 9 a > d <0 => 0<cM<\C g \ => \C x \ = \C g \-CM<\C B \ ^>1. (E2) 



APPENDIX F: COMBINED HEAT CAPACITY Cx + C^ } AS A FUNCTION OF R g 

This appendix shows a behaviour of combined heat capacity Cx + Cy* 1 as a function of black hole radius R g . We 
make use of the calculations done in appendix [Cj By equation (|C3[) and definitions of Cx and Cy given in equations 
(fTTj) , Cx + Cy^ is expressed as 

Cx(R g ) + Cy"\R g ) = -2irR 2 h x 2 + ^ [f(x)+qx 3 (x+ Vl-a; 2 )] , (Fl) 

where g = (ATrRhTh) 4 , x = R g /Rh, < x < 1 by definition, and /(x) is given by equation (|C4|) . Then, using the 
limit values (IC5D . we find 

x) as x — > 

, a' a' ■ (F2) 

2-kRi + + t-^t Q at x = 1 




8tt 2 8tt 2 



The first two terms in Cx + Cy ' 1 at x = 1 is negative in the framework of NPT model as discussed in equations (|C5[) . 
Therefore, if Th is sufficiently small, then q can be small enough so that Cx + Cy^ at x = 1 becomes negative. 

The differential of Cx + Cy^ is complicated and not suitable for analytical utility. Instead of analytical discussion, 
we show some numerical examples of Cx + Cy* 1 in figure [51 It is recognised with this figure that, for a sufficiently 
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FIG. 8: Numerical plots of Cx + C\?> as a function of R g with N = 10. Upper three panels are for Rh = 2 and 7), = 0.5, 
0.27 and 0.1 from left to right. Lower three panels are for Rh = 10 and Th = 0.05, 0.0013 and 0.001 from left to right. It has 
been found for every value of Rh > 1 as far as the author checked that the same behaviour is observed with decreasing Th, and 
Cx + C ( y j) < can hold for a sufficiently small Th- 



small Th, the condition (fT2")l . Cx + Cy < 0, which guarantees the validity of NPT model can be satisfied for a 
certain range of R g . Although figure [S] shows only some examples, it has been found as far as the author checked that 

the same behaviour is observed for every value of R h > 1. That is, C x + C^ ] < can hold for a sufficiently small 
Th- Here note that, during a semi-classical and quasi-equilibrium stage of black hole evaporation, it is reasonable to 
require Th < 1 (= Planck temperature ). Hence we assume throughout the paper that Th is small enough so that the 

equation Cx(R g ) + Cy (R g ) = has one, two or three solutions of R g . 

Finally in this appendix, we discuss what happens along a black hole evaporation process for the case that the 
equation Cx{R g ) + Cy\R g ) = has two or three solutions of R g . For the first consider the case of three solutions, 
and denotes these solutions by R g Q , R^ and R g Q in an increasing order {R g o < RIq < R g o ) ■ When the evaporation 

starts with initial black hole radius R g {0) which is in the range R g l Q < R g (Q) < Rh (see for example upper-centre 
panel in figure [5]), then the evaporation process can be treated in the framework of NPT model until the radius R g {t) 
decreases to RgQ- According to the discussion in section Ull Bl (the paragraph one before the last one), the evaporation 
process can not be treated in the framework of NPT model after the black hole radius decreases to R g o, because 

the validity condition (fT2"|) of NPT model is violated in the range R g m ^ < R g < R g o- However it is reasonable to 
expect that, even if the NPT model is not applicable to analyse the evaporation process while the black hole radius 
decreases from Rjp to Rg™\ the NPT model becomes valid again after the black hole radius reaches R g ™\ That is, 
the evaporation process can be analysed in the framework of NPT model in the range R g S Q < R g < R g ^ - 

Next consider the case that the equation Cx{R g )+Cy {R g ) = has two solutions. When the black hole evaporation 

starts with initial radius which is larger than the larger solution of Cx{R g ) + Cy {R g ) = (see for example lower- 
centre panel in figure [5]), then the similar discussion as given in previous paragraph can be applied. That is, although 
the black hole evaporation process can not be treated in the framework of NPT model until the black hole radius 
becomes smaller than the larger solution of Cx{R g ) + Cy\R g ) = 0, it becomes possible that the evaporation process 
can be analysed by the NPT model until the radius decreases to the smaller solution of Cx{R g ) + Cy \R g ) = 0. 
From the above we find that, for both of the cases that the equation Cx{R g ) + Cy\R g ) = has two and three 

solutions, it is the lowest solution of the equation C x {R g ) + Cy{R g ) = at which the evaporation process comes 
not to be treated in the framework of NPT model and the evaporation process proceeds to a highly nonequilibrium 
dynamical stage. (According to section Ull C[ the highly nonequilibrium dynamical stage is the abrupt catastrophic 
evaporation at semi-classical level for the case Rh > O(10 3 ).) 
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APPENDIX G: ON THE FINAL STAGE OF EVAPORATION: EXCEPT FOR THE NPT MODEL 



This appendix shows that the abrupt catastrophic evaporation at semi-classical level does not occur for the equi- 
librium model used in section HT1 and for the model that a black hole evaporates in an empty space. 

For the first consider the equilibrium model used in section [TXJ Since there is no nonequilibrium region between 
black hole and heat bath, the energy transport equations are given by 



dEg 

dt 



' rp4 



A„ 



dE h 
dt 



a'(T £ 



A„ 



(Gl) 



where A g = 4irR 2 is the area of black hole surface. 
Then one of equations of states (P), E g = R g /2, gives 



The first equation of dE g /dt is used in the following analysis. 



dRg 


= 2 


dE g 


dt 




dt 



2a'{T^-n)A g . 



(G2) 



Here note that, since the equilibrium model is based on equations of states (JlJ, the quasi-equilibrium assumption is 
also necessary for this model and v < 1 is required. It is obvious that v = 1 occurs for non-zero radius R g > 0, since 
T g A g cx 1/-Rg — > oo as R g — ► 0. The inequality v < 1 corresponds to the validity condition (TT3"]) of NPT model. On 
the other hand, for the heat capacity C g of black hole, it is obvious that C g = R g = 0. That is, there is no 



counterpart in the equilibrium model which corresponds to the validity condition (|12p of NPT model. Hence we can 
recognise that, if v = 1 occurs for a black hole of semi-classical size R g > 1, it is reasonable to conclude that the 
abrupt catastrophic evaporation at semi-classical level occurs. And if v = 1 does not occur for a semi-classical black 
hole, then the abrupt catastrophic evaporation does not occur in the framework of equilibrium model. 
The validity condition v < 1 is rewritten as follows, 



< 2tt Tt 



rpl 



R^ 



l 



8tt 2 /3 



1 



1 + (2T, 2 //?) 



(G3) 



where (3 = 2tt/<t' = 240/irN and R g = l/47rT g is used. Here equation (J3J| gives /? = O(10). Further note that T g 
holds generally for any evaporation process and that T g < 1 holds due to the quasi-equilibrium assumption. Then we 
find Th < 1 and (2T%/(3) 2 <C 1. Hence we can approximate inequality (|G3[) to the following form 



R 2 g> 



l 



8tt 2 P 



Rg £ 0(0.01) . 



(G4) 



where (3 = O(10) is used. This denotes that the fast evaporation v = 1 occurs at R g ~ 0.01. Hence, as discussed 
in previous paragraph, it is reasonable to conclude that the abrupt catastrophic evaporation does not occur in the 
framework of equilibrium model used in section HT1 

Next consider a black hole evaporation in an empty space. Ignoring gravitational interactions between black hole 
and radiation fields (curvature scattering, lens effects and so on), the Stefan-Boltzmann law gives (see the energy 
emission rate J empty defined by equation (|17p) 



Then one of equations of states (Q]), E g — R g /2, gives 



dt - aT ° Ag 



dRg 


= 2 


dEg 


dt 




dt 



= 2a'T^A g . 



Hence, following a similar calculation given in previous paragraph with setting Th 

1 



0, we obtain 



v < 1 



0{lQ- 6 ) 



R g > 0(0.01) . 



(G5) 



(G6) 



(G7) 



9 8tt 2 /3 

This denotes that the fast evaporation v = 1 occurs at R g ~ 0.01, and that it is reasonable to conclude that the abrupt 
catastrophic evaporation does not occur for a black hole evaporation in an empty space with ignoring gravitational 
interactions between black hole and radiation fields. 
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